DUDLEY  KNOX  LIB;     , 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY  CA    93943-5101 


Approved  for  public  release;  distribution  is  unlimited. 

Active  Damping  Control  of  a  Flexible  Space  Structure  using 
Piezoelectric  Sensors  and  Actuators 

by 

Scott  M.  Newman 

Lieutenant,  United  States  Navy 

B.S.  Aerospace  Engineering,  United  States  Naval  Academy, 

1984 

Submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  ASTRONAUTICAL  ENGINEERING 

from  the 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 
Unclassified 


lb  RESTRICTIVE  MARKINGS 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION/DOWNGRADING  SCHEDULE 


3  DISTRIBUTION/AVAILABILITY  OF  REPORT 
Approved  for  public  release;  distribution  is  unlimited. 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 
Naval  Postgraduate  School 


6b  OFFICE  SYMBOL 
(If  applicable) 
55 


7a  NAME  OF  MONITORING  ORGANIZATION 
Naval  Postgraduate  School 


6c.  ADDRESS  (City.  State,  and  ZIP  Code) 
Monterey,  CA  93943-5000 


7b  ADDRESS  (City.  State,  and  ZIP  Code) 
Monterey,  CA  93943-5000 


8a.  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 


8b  OFFICE  SYMBOL 
(If  applicable) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


8c  ADDRESS  (City,  State,  and  ZIP  Code) 


10.  SOURCE  OF  FUNDING  NUMBERS 


Program  Element  No 


Project  No 


Task  No 


Work  Unit  Accession 
Number 


1 1 .  TITLE  (Include  Security  Classification) 

Active  Damping  Control  of  a  flexible  Space  Structure  Using  Piezoelectric  Sensors  and  Actuators 


1  2   PERSONAL  AUTHOR(S)    Newman,  Scott  M. 


13a.  TYPE  OF  REPORT 
Master's  Thesis 


13b.  TIME  COVERED 
From  To 


14  DATE  OF  REPORT  (year,  month,  day) 
December  1992 


15.  PAGE  COUNT 
74 


16.  SUPPLEMENTARY  NOTATION 

The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy  or  position  of  the  Department  of  Defense  or  the  U.S. 

Government. 


17.  COSATI  CODES 


FIELD 


GROUP 


SUBGROUP 


18  SUBJECT  TERMS  (continue  on  reverse  if  necessary  and  identify  by  block  number) 

Active  Damping  Control,  Piezoelectric  Sensors  and  Actuators,  Positive  Position  Feedback, 
Strain  Plate  Feedback,  Flexible  Spacecraft  Simulator  (FSS) 


1 9.  ABSTRACT  (continue  on  reverse  if  necessary  and  identify  by  block  number) 

This  thesis  details  the  experimental  analysis  of  an  active  damping  control  technique  applied  to  the  Naval  Postgraduate  School's  Flexible 
Spacecraft  Simulator  using  piezoceramic  sensors  and  actuators.  The  mass  property  of  the  flexible  arm  is  varied  to  study  the  frequency  effects  on 
the  Positive  Position  Feedback  (PPF)  algorithm.  Multi-modal  dynamic  response  is  analytically  studied  using  a  finite-element  model  of  a  canti- 
levered  beam  while  under  the  influence  of  three  different  control  laws:  a  basic  law  derived  from  the  Lyapunov  Stability  Theorem,  PPF  and 
Strain  Rate  Feedback  (SRF).  The  advantages  and  disadvantages  of  using  PPF  and  SRF  for  active  damping  control  are  discussed. 


20.  DISTRIBUTION/AVAILABILITY  OF  ABSTRACT 

Q  UNCLASSIFIED/UNLIMITED  Q   SAME  AS  REPORT  ]    DTIC  USERS 


21.  ABSTRACT  SECURITY  CLASSIFICATION 
Unclassified 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 
Agrawal,  Brij  N. 


22b  TELEPHONE  (Include  Area  code) 
(408)656-3338 


22c  OFFICE  SYMBOL 
AA/Ag 


DD  FORM  1473,  84  MAR 


83  APR  edition  may  be  used  until  exhausted 
All  other  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

Unclassified 


ABSTRACT 

This  thesis  details  the  experimental  analysis  of  an  active 
damping  control  technique  applied  to  the  Naval  Postgraduate 
School's  Flexible  Spacecraft  Simulator  using  piezoceramic  sen- 
sors and  actuators.  The  mass  property  of  the  flexible  arm  is 
varied  to  study  the  frequency  effects  on  the  Positive  Position 
Feedback  (PPF)  algorithm.  Multi-modal  dynamic  response  is 
analytically  studied  using  a  finite- element  model  of  a  canti- 
levered  beam  while  under  the  influence  of  three  different 
control  laws:  a  basic  law  derived  from  the  Lyapunov  Stability 
Theorem,  PPF  and  Strain  Rate  Feedback  (SRF) .  The  advantages 
and  disadvantages  of  using  PPF  and  SRF  for  active  damping 
control  are  discussed. 


111 


TABLE  OF  CONTENTS 

I.  INTRODUCTION   1 

A.  BACKGROUND 1 

B.  FOCUS  OF  THESIS 4 

II.  THEORETICAL  ANALYSIS  5 

A.  PIEZOELECTRIC  THEORY   5 

B.  FINITE  ELEMENT  MODEL   8 

1.  Equations  of  Motion   8 

2 .  Actuator  Moment 13 

3 .  Modal  Coordinate  Form 15 

C.  COMPUTER  SIMULATIONS 18 

1.  Lyapunov  Stability-Derived  Control  Law  ...  19 

2.  Positive  Position  Feedback  (PPF)  24 

3.  Strain  Rate  Feedback  (SRF) 34 

4 .  Comparison  of  PPF  and  SRF 44 

III.  EXPERIMENTAL  ANALYSIS 48 

A.  PHYSICAL  SET-UP 48 

B.  PROCEDURE 50 

C.  RESULTS 52 

IV.  CONCLUSIONS 58 


IV 


>U=YK,N 


APPENDIX 60 


REFERENCES 64 


BIBLIOGRAPHY   65 


INITIAL  DISTRIBUTION  LIST  67 


ACKNOWLEDGMENT 

The  author  would  like  to  express  his  gratitude  and  appre- 
ciation to  Dr.  Hyochoong  Bang  for  the  endless  hours  of  his  tu- 
toring, insight  and  encouragement.  The  progress  of  this  thesis 
would  not  have  been  achieved  without  his  assistance  and  exper- 
tise. An  equal  amount  of  thanks  goes  to  Dr.  Brij  N.  Agrawal 
for  providing  the  resources  used  with  this  project  and  for 
guiding  the  experiment.  His  assistance  with  defining  the  focus 
of  the  thesis  and  his  encouragement  during  difficult  times 
were  invaluable. 


VI 


I .   INTRODUCTION 

A.   BACKGROUND 

The  evolution  of  spacecraft  into  larger  and  more  complex 
structures,  such  as  Space  Station  Freedom,  has  greatly  expand- 
ed the  problem  of  structural  dynamics.  The  need  to  limit 
launch  costs  drives  the  engineer  to  a  low  mass  design  which 
results  in  these  large  structures  being  extremely  flexible 
with  very  low  frequency  fundamental  vibrational  modes.  These 
structural  characteristics  pose  some  difficult  control  prob- 
lems which  are  the  focus  of  the  emerging  Control -Structures 
Interaction  (CSIJ  field. 

There  are  a  variety  of  scenarios,  such  as  slewing  and 
pointing  maneuvers,  docking  with  other  spacecraft  and  inter- 
actions with  environmental  forces  that  will  produce  unwanted 
structural  vibrations.  There  may  be  requirements  for  sensor  or 
antenna  pointing  or  attitude  control  which  can  not  be  met  with 
only  the  inherent  structural  damping  in  effect.  This  describes 
the  need  to  devise  techniques  which  actively  control  the  dy- 
namic response  of  flexible  structures.  Any  technique,  in  sim- 
plistic terms,  is  composed  of  two  separate,  but  interrelated 
operations.  First,  the  motion  of  the  structure  must  be  accur- 
ately measured  by  some  type  of  sensor.  Then,  this  information 


must  be  fed  back  via  an  effective  control  algorithm  to  an  ac- 
tuator on  the  structure  that  will  counteract  the  vibrations. 

When  trying  to  measure  the  dynamics  of  a  structure,  two 
possible  quantities  that  can  be  measured  are  the  acceleration 
of  a  point  or  points  on  the  structural  member  and  the  strain 
within  the  member.  However,  with  the  low  frequency  vibrations 
characteristic  of  the  structures  we  are  interested  in,  precise 
acceleration  measurements  are  difficult.  This  leaves  strain  as 
the  quantity  best  suited  for  defining  low  frequency  vibra- 
tions . 

There  are  a  variety  of  sensors  available  for  this  applica- 
tion such  as  fiber  optic  sensors,  strain  gages  and  piezoelec- 
tric ceramics.  Fiber  optic  sensors  can  be  used  to  make 
measurements  with  a  variety  of  methods,  but  are  not  yet  prac- 
tical for  complexity  and  technology  problems.  Strain  gages  are 
limited  by  a  relatively  low  sensitivity  and  are  highly  subject 
to  noise.  The  one  type  that  has  received  much  study  and  is 
seeing  increased  usage  is  the  piezoelectric  ceramic.  Reference 
1  highlights  its  advantages  as  having: 

•  high  strain  sensitivity 

•  low  noise 

•  low  to  moderate  temperature  sensitivity 

•  easy  implementation 

The  desired  characteristics  of  a  vibration  suppression  ac- 
tuator are  low  mass,  infinite  bandwidth,  electrically  powered 
and  an  internal  force  producer  [Ref .  2] .  A  low  mass 


requirement  is  necessary  so  as  not  to  significantly  alter  the 
physical  characteristics  of  the  structure  as  well  as  the 
overall  launch  mass.  An  ideal  infinite  bandwidth  would  theor- 
etically allow  control  of  all  vibrational  modes  of  the  struc- 
ture. Electric  power  reduces  the  complexity  and  mass  inherent 
in  a  mechanical  system.  An  actuator  that  produces  internal 
forces  only  would  greatly  simplify  the  control  algorithm  for 
the  entire  structure. 

A  number  of  materials  have  been  developed  for  this  purpose 
such  as  piezoelectric  ceramic  wafers  (PZT) ,  electrorestrictive 
ceramic  wafers  (PMN) ,  piezoelectric  polymer_  film  (PVF2)  and 
shape  memory  metal  (nitinol)  wires.  PZT  and  PMN  can  produce 
similar  strain  levels,  but  PMN  is  non- linear  and  very  temper- 
ature dependent.  Although  nitinol  can  produce  strains  3  00 
times  greater  than  PZT,  it  is  difficult  to  incorporate  and  not 
practical  for  most  applications.  When  considering  these  trade- 
offs, PZT  becomes  the  best  choice  for  an  actuator.  Its  suita- 
bility is  summarized  in  the  following  characteristics  [Ref .1] : 

•  high  stiffness 

•  sufficient  stress  to  control  vibrations 

•  good  linearity 

•  temperature  insensitivity 

•  easy  to  implement 

•  low  power  consumption 


To  further  investigate  active  damping  control  techniques, 
this  technology  has  been  incorporated  into  the  Naval  Postgrad- 
uate School's  Flexible  Spacecraft  Simulator  (FSS) . 

B.   FOCUS  OF  THESIS 

The  effectiveness  and  utility  of  using  piezoelectric  sen- 
sors and  actuators  to  control  vibrations  will  be  analytically 
demonstrated  via  computer  simulations.  A  model  will  be  devel- 
oped for  a  simple  cantilevered  beam  using  the  finite  element 
method.  The  dynamic  response  of  this  system  will  be  analyzed 
while  under  the  influence  of  three  different  control  laws: 

•  a  basic  Lyapunov  stability- derived  control  law 

•  Positive  Position  Feedback  (PPF) 

•  Strain  Rate  Feedback  (SRF) 

The  Positive  Position  Feedback  algorithm  will  also  be  in- 
vestigated experimentally.  This  is  the  particular  control  law 
which  has  been  built  into  the  Flexible  Spacecraft  Simulator. 
The  mass  property  of  the  flexible  arm  can  be  adjusted  to  vary 
the  natural  frequency  of  the  first  vibrational  mode.  As  will 
be  shown,  the  effect  of  the  piezoelectric  actuator  on  the  dy- 
namic response  of  the  beam  will  change  as  the  natural  frequen- 
cy varies  from  a  particular  optimum  value.  This  will  highlight 
the  advantages  and  disadvatages  of  this  technique. 


II.   THEORETICAL  ANALYSIS 

A.   PIEZOELECTRIC  THEORY 

Piezoelectric  ceramics  and  crystals  derive  their  utility 
in  vibration  control  from  the  ability  to  convert  electrical 
energy  into  mechanical  energy  and  vice  versa.  When  a  force  is 
applied  to  the  material,  the  crystalline  structure  creates  a 
voltage  proportional  to  the  force.  Likewise,  an  electric  po- 
tential applied  across  the  material  will  alter  its  shape 
dimensional ly. 

Man-made  ceramics  have  an  advantage  over  natural  crystals 
in  that  they  can  be  customized  for  each  application.  These 
ceramics  can  be  manufactured  in  almost  any  size  and  shape  and 
have  their  piezoelectric  properties  oriented  in  any  specified 
direction.  These  properties  are  induced  into  the  material  dur- 
ing a  process  called  "poling".  While  heated  to  a  temperature 
above  its  Curie  point,  an  electric  field  is  applied  across  the 
material.  This  has  the  effect  of  realigning  the  crystalline 
structure  such  that  it  is  elongated  in  a  direction  parallel  to 
the  electric  field,  the  "poling"  axis,  and  foreshortened  in  a 
direction  perpendicular  to  the  "poling"  axis.  Henceforth,  as 
shown  in  Figure  1,  any  applied  voltage  of  the  same  polarity  as 
the  "poling"  voltage  will  elongate  the  ceramic  along  the 
"poling"  axis  and  a  voltage  of  opposite  polarity  will  elongate 
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Figure  1  Piezoelectric  actuator  mode 


it  in  a  direction  perpendicular  to  the  "poling"  axis.  In  the 
reverse  operation,  a  parallel  compressive  force  or  a  perpen- 
dicular tensile  force  in  relation  to  the  "poling"  axis  will 
produce  a  voltage  of  the  same  polarity  and  vice  versa.  (See 
Figure  2 . )  [Ref .  3] 

The  voltage  (Vs)  produced  by  a  sensor  is  given  by  the 
relationship: 
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where  E  is  the  modulus  of  elasticity  of  the  ceramic,  d3,  is  the 
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Figure  2  Piezoelectric  sensor  mode 


lateral  strain  coefficient,  D  is  the  dielectric  constant,  t  is 
the  thickness  and  el  and  e2  are  the  lateral  strains  [Ref .4] . 
This  lays  the  groundwork  for  measuring  the  vibration  of  a 
structural  member  by  converting  the  strain  on  the  ceramic  pro- 
duced by  the  flexural  vibrations  into  a  voltage  which  can  be 
much  more  easily  manipulated  and  processed. 

The  effect  of  an  actuator  is  described  by  the  equation  for 
the  bending  moment  (mb(x,t))  at  a  cross- section  of  the  struc- 
tural member: 


mb(x,t)    =EIy"(x,t)    -  kdvAt)  (2) 


where  EI  is  the  flexural  rigidity  of  the  member,  y"(xft)  is 
the  second  order  partial  derivative  of  the  transverse  dis- 
placement with  respect  to  x  (the  longitudinal  coordinate) ,  kd 
is  a  constant  dependent  on  the  ceramic  piezoelectric  constant, 
dielectric  constant,  width  and  the  member  thickness  and  vd(t) 
is  the  voltage  applied  to  the  ceramic  [Ref .  5]  .  This  shows  how 
an  applied  voltage  can  produce  a  strain  which,  in  effect, 
places  concentrated  moments  (a  couple)  at  the  endpoints  of  the 
element  on  which  the  piezoceramic  material  is  located.  The 
direction  of  the  moment  needed  to  counteract  a  vibration  is 
changed  by  changing  the  polarity  of  the  applied  voltage. 

B.   FINITE  ELEMENT  MODEL 

In  order  to  effectively  analyze  the  dynamic  response  of  a 
system  and  the  effect  of  an  applied  control  law,  an  accurate 
and  workable  model  of  the  system  must  be  developed.  The  finite 
element  method  is  used  to  discretize  the  equations  of  motion 
from  unwieldy  partial  differential  equations  into  second-order 
differential  equations  of  generalized  coordinates. 

1.   Equations  of  Motion 

This  analysis  will  consider  a  cantilevered  beam  as  the 
system.  From  Bernoulli-Euler  beam  theory,  the  kinetic  energy 
(T) ,  potential  energy  (U) ,  and  work  (W)  done  by  an  applied 


moment  from  a  piezoceramic  actuator  (M(x) )  on  each  element 
are: 
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where  x  is  the  coordinate  in  the  longitudinal  direction,  y  is 
the  coordinate  in  the  transverse  direction,  Ve  is  the  trans- 
verse velocity,  E  is  the  modulus  of  elasticity,  I  is  the  mo- 
ment of  inertia  about  the  bending  axis,  p  is  the  mass  per  unit 
length,  Xj  is  the  coordinate  of  the  left  end  of  the  element  and 
h  is  the  length  of  the  element.  A  typical  finite  element  is 
shown  in  Figure  3  with  the  associated  moment  (M)  and  shear  (V) 
forces . 

From  the  Extended  Hamilton's  Principle,  the  following 
relationship  is  used  to  formulate  the  equations  of  motion: 
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8  (T-U+W)    dt  =   0  (6) 


Taking  a  variation  of  each  term  and  assuming  a  constant  moment 


over  the  length  of  the  element  yields: 
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Figure  3  Finite  Element  Schematic 

Integrating  the  first  term  by  parts  gives  the  following  varia^ 
tional  equation: 

dt  =   0 
(8) 


Jt,    Jx,   K      dx2{  dx2  Jx,   ax2 
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At  this  point,  substitutions  are  made  to  discretize  the  motion 
into  the  form: 


y{x,t)    =  J2  NiM   9Ti(t) 


(9) 


The  N,  terms  are  third- order  polynomials,  chosen  to  satisfy  the 
boundary  conditions  at  both  ends  of  the  element,  which  define 
the  mode  shapes  and  are  functions  of  x  only  [Ref .  6] : 
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The  qj  terms  are  the  system  generalized  coodinates  which  are 
functions  of  time  only.  The  quantities  they  represent  are: 


Transverse  Deflection  at  Left  End  of  Element 
Rotation  at  Left  End  of  Element  (11) 

Transverse  Deflection  at  Right  End  of  Element 
Rotation  at  Right  End  of  Element 


Taking  the  appropriate  derivatives  of  Equation  (9)  and  substi- 
tuting back  into  Equation  (8)  and  taking  a  summation  of  the 
ith  terms  yields: 
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(12) 


Since  this  relationship  must  hold  true  over  any  time  period 
and  for  any  arbitrary  variation  dq7,  the  integrand  must  equal 
zero.  Now,  the  familiar  second- order  equation  of  a  vibratory 
system  becomes  evident  in  the  form: 


M^q^t)    +K1Jqi(t)    =AG^(x)£ 


x^h 


,  J  =  1..J2 


(13) 


The  mass  (M)  and  stiffness  (K)  matrices  are  computed  from: 


M£j  =  p  |Xi+    2^  (x)  Nj  (x)    dx 
K„  =  EI  fXi+h  N"(x)N"(x)    dx 


(14) 


Performing  the  integrations  results  in  a  system  of  four  equa' 
tions  for  each  element: 
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The  complete  mass  and  stiffness  matrices  for  the  entire  system 
are  obtained  from  a  simple  combination  of  the  elemental  matri- 
ces. Since  the  q$  and  q4  coordinates  of  one  element  are  the 
same  as  the  q,  and  q2  coordinates  of  the  next  element,  the  glo- 
bal matrix  is  formed  by  overlapping  (adding)  the  upper  left 
2x2  submatrix  of  an  element  with  the  lower  right  2x2  submatrix 
of  the  previous  element.  Therefore,  a  n- element  system  has  ma- 
trices with  dimensions  of  (nx2)+2  by  (nx2)+2.  However,  for  a 
cantilevered  beam,  the  qj  and  q2  coordinates  of  the  first  ele- 
ment are  always  zero  reducing  the  matrices  to  nx2  by  nx2 . 
2 .   Actuator  Moment 

The  control  moment  produced  by  an  actuator  on  an  ele- 
ment is  derived  using  the  Lyapunov  Stability  Theorem.  The  to- 
tal energy  of  the  beam: 


2  Jx, 
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dx 


(16) 


is  chosen  as  the  Lyapunov  function.  This  function  meets  the 
requirements  of  being  positive  definite,  continuously 
dif ferentiable  and  equal  to  zero  at  the  equilibrium  state.  The 
only  remaining  condition  needed  to  ensure  stability  is  that 
the  time  derivative  of  U  must  be  less  than  or  equal  to  zero. 
[Ref.  7] 
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The  time  derivative  of  the  energy  function  is: 


u  =  f 


oyy  +  ei — —  —  — — 
Pyy  dx*dt[dx>) 


dx 


(17) 


Integrating  this  equation  using  the  following  relationships 
for  a  vibrating  beam,  bending  moment  (M)  and  shear  (V) : 


py  +  £J-^-Z  =  0 


(18) 


M  =  EI 


&y 
dx2 


(19) 


V  =    *"  -  EI*? 
dx  dx2 


(20) 


while  assuming  no  shear  in  the  beam  yields 


U  =  M  (^   -  ^ 


(21) 


The  quantity  in  parentheses  can  be  measured  by  a  specially 
designed  electric  circuit  which  amplifies  the  current  devel- 
oped over  the  sensor.  To  ensure  that  the  derivative  of  the 
Lyapunov  function  is  always  less  than  zero,  M  must  be  of  the 
form: 
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M=  -k  (yii+h  -  yij 


(22) 


M  is  directly  proportional  to  the  applied  voltage  to  the  actu 
ator  such  that  k  is  a  positive  constant  resulting  in: 


v  =  k  (y£1+i  -  yi,)2  ±  o 


(23) 


In  terms  of  the  generalized  coordinates,  the  moment  is  ex- 
pressed as: 


M  =  -k  (g4  -  q2) 


(24) 


The  system  of  equations  for  an  element  with  an  actuator  at 
tached  now  becomes: 


Mij%i  +  KijQi  = 


0  0  0  0  1  Si- 

o  -k  o  k   $2 

0     0  0  0     (j3 

0  k  0  -Jfej  q 


(25) 


3 .   Modal  Coordinate  Form 

The  equations  of  motion  can  be  decoupled  and  more  eas- 
ily analyzed  by  transforming  the  generalized  coordinates  into 
modal  coordinates.  This  transformation  is  of  the  form: 


q  =  *l 


(26) 
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where  $  is  the  weighted  modal  matrix  and  £  is  a  vector  of 
modal  coordinates . 

First,  the  eigenvalues  and  eigenvectors  must  be  calcu- 
lated from  the  freely  vibrating  undamped  system: 

Afg  +  Kq  =  0  (27) 

For  a  n-dimensional  system,  there  are  n  eigenvalues  (X,-,  i= 
l...n)  and  n  eigenvectors  (#,,  1=1... n)  which  satisfy  the 
equation: 

XjMbi   =  K$i  (28) 

The  eigenvectors,  also  known  as  the  normal  modes,  are  orthog- 
onal with  respect  to  the  mass  and  stiffness  matrices.  This 
property  allows  the  matrices  to  be  diagonalized  using  the  re- 
lationships : 

$1*®!  =  «i  (29) 

where  the  diagonal  elements  M^  and  K-,  (i=l...n),  are  known  as 
the  generalized  mass  and  stiffness  elements  for  each  mode. 

The  weighted  modal  matrix  (#)  is  obtained  by  dividing 
the  ith  column  of  the  modal  matrix  by  the  square  root  of  the 
ith  generalized  mass: 
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+ 


iJ 


i=l..j3  ,  j=l...n 


(30) 


Substituting  the  weighted  modal  matrix  into  Equation  (29)  pro- 
duces the  diagonalizations : 


<&TM$   =  I 


(31) 


where  I  is  the  identity  matrix  and  the  diagonal  elements  of  K* 
are  the  eigenvalues  and  squares  of  the  modal  frequencies. 

If  the  system  is  assumed  to  have  viscous  damping  and 
excited  by  an  arbitrary  function  Fu,  the  equations  of  motion 
take  the  form: 


Mq  +  Cq  +  Kq   =  Fu 


(32) 


If  the  damping  matrix  is  proportional  to  either  M  or  K;  it  is 
also  diagonalized  by  the  weighted  modal  matrix.  Applying  the 
transformation  of  Equation  (26)  to  Equation  (32)  and  premul- 
tiplying  by  *  results  in  n-uncoupled  equations  of  the  form: 


I 


2C« 


i 


<*>' 


I  =  ®TFu 


(33) 
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C.   COMPUTER  SIMULATIONS 

To  demonstrate  the  utility  and  stability  of  the  specified 
control  laws,  the  finite  element  model  has  been  programmed 
using  MATLAB.  The  basic  approach  is  to  put  the  modal  equations 
into  a  state- space  form  and  use  MATLAB' s  linear  simulator 
function  to  produce  a  time  response.  The  code  for  the  PPF 
simulation  is  in  the  Appendix.  The  code  for  the  other  control 
laws  is  of  a  similar  form  except  for  the  specific  control  law 
formulation. 

The  following  values  for  the  physical  parameters  of  the 
cantilevered  beam  were  selected: 

•  number  of  elements:  n  =  8 

•  element  length:  h  =  .1  m 

•  element  density:  p    =    .5  kg/m3 

•  product  of  modulus  of  elasticity  and     EI  =  .9  Nm2 
moment  of  inertia 

Using  these  values,  the  system  mass  and  stiffness  matrices 

were  formulated  and  used  to  solve  the  eigenvalue  problem.  That 

solution  provided  the  modal  frequencies  of  the  beam.  The  first 

three  are: 

•  coj  =  7.37  rad/s  =  1.17  Hz 

•  co2  =  46.19  rad/s  =  7.35  Hz 

•  co3  -  129.42  rad/s  =  20.6  Hz 

The  inherent  structural  damping  was  introduced  by  arbitrarily 
damping  ratios  (£",.)  such  that  the  damping  increased  from  the 
lowest  value  for  the  first  mode  to  the  highest  value  for  the 
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nth  mode.  In  general,  the  high  frequency  modes  are  damped  out 
more  quickly  then  the  lower  frequency  modes  due  to  the  inher- 
ent natural  frequencies. 

For  all  the  simulations,  the  beam  was  deflected  and  re- 
leased from  rest  such  that  the  first  three  modes  were  excited. 
Simple  geometric  calculations  were  performed  to  determine  the 
required  initial  values  for  each  of  the  generalized  coordi- 
nates. The  open- loop  response  of  the  vibrating  beam  is  shown 
in  Figure  4 . 

From  this  point  on,  the  simulations  are  customized  to  the 
specific  control  law  utilized.  The  simulation  development  and 
results  are  detailed  in  the  following  sections  for  each  of  the 
three  chosen  control  laws. 

1.   Lyapunov  Stability-Derived  Control  Law 

As  shown  in  Chapter  II. B. 2,  the  equations  of  motion 
for  each  element  are  of  the  form: 

Mifii  +  Ku<Ii  =  Fk*Zi  (34) 

where  k*  is  a  1  x  4  vector  containing  the  proportionality 
factor  k: 

k*  =   [  0  k  0  -Jfc  ]  <35) 

F  is  a  4  x  1  vector  dependent  on  actuator  placement.  If  an 
actuator  is  located  on  an  element, 
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Figure  4  Open- Loop  Response 


F  =  [  0  -1  0  1  ] 


(36) 


Otherwise,  it  is  a  null  vector.  For  the  entire  system  of  n 
elements,  the  individual  4x4  Fk*  matrices  are  combined  to 
form  the  n  x  n  matrix  F*.  In  modal  form,  the  system  is 
represented  by: 


I    +    (D  -  F)l   +  K*Z   -  0 


(37) 


where : 
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2d"! 


»2 


D  = 


2C2« 


2W2 


f  =  &tf  *<&   jr 


to- 


2C  w 


to, 


The  system  is  further  manipulated  into  state  space  form: 


X  =  AX  +  Bu 


Y   =  CX  +  Du 


(38) 


where 


*-[**]' 


A  = 


0     I 
K*   F  -  D 


C  is  the  identity  matrix  and  B  and  D  are  null  matrices. 

Four  different  simulations  were  run  with  varying  actu- 
ator locations  as  follows: 

•  Case  1 .     Actuator  at  base  (element  1) 

•  Case  2 .     Actuator  in  middle  (element  5) 

•  Case  3 .     Actuator  at  tip  (element  8) 

•  Case  4 .  Combination  of  Cases  1-3  (elements  1,5  and  8) 
The  results  are  shown  in  Figures  5-8.  As  expected,  the  actu- 
ator at  the  base  has  the  most  significant  effect  on  damping 
vibrations  since  it  can  apply  the  largest  moment.  Although 
actuators  in  the  middle  and  at  the  tip  do  provide  significant 
damping  when  compared  to  the  open- loop  response,  their  contri- 
bution is  minimal  in  the  combined  system.  A  comparison  of  Case 
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1  and  Case  4  shows  how  much  the  base  actuator  dominates.  This 
fact  highlights  the  advantage  of  using  a  single  actuator  when 
considering  cost  and  implementation  issues. 
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Figure  5  Lyapunov  /  Actuator  at  base 

The  proportionality  factor  k  is  a  function  of  the 
physical  constants  of  the  piezoelectric  which  relate  the  ap- 
plied voltage  to  the  moment  created  and  a  gain  factor  applied 
to  the  feedback  circuit.  For  simplicity,  it  is  assumed  to  be 
unity  for  these  simulations. 

Varying  k  has  a  significant  effect  on  the  response. 
Theoretically,  the  system  is  stable  for  any  positive  value. 
However,  larger  values  are,  in  effect,  placing  a  moment  with 
a  magnitude  larger  than  required  which  tends  to  "overcontrol " 
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Figure  6  Lyapunov  /  Actuator  in  middle 
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Figure  7  Lyapunov  /  Actuator  at  tip 
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Figure  8  Lyapunov  /  All  three  actuators 

the  beam.  This  is  seen  in  the  response  as  the  actuators,  al- 
though stabilizing,  drive  the  beam  to  excessive  amplitudes. 
Small  values  of  k  equate  to  moments  of  insufficient  magnitude 
which  decrease  the  damping  ratio.  Optimizing  the  system  re- 
quires a  tradeoff  between  the  desired  settling  time  and  a  max- 
imum allowable  amplitude.  A  trial  and  error  procedure  is 
required  to  select  the  best  value  and  customize  the  control  to 
the  application. 

2.   Positive  Position  Feedback  (PPF) 

This  control  algorithm  is  defined  strictly  in  the 
modal  form.  For  a  scalar  case,  the  equations  of  motion  are 
represented  by: 
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£    +  2C&>£    +  <o2£   =  Gu>2T\ 


(39a) 


il    +  2Cc<octi    +  o)Jii   =  (0^ 


(39b) 


where  £  is  the  structure  modal  coordinate,  r\  is  the  compensa- 
tor modal  coordinate,  $  and  j*c  are  the  structure  and  compensa- 
tor damping  ratios,  to  and  toc  are  the  structure  and  compensator 
natural  frequencies  and  G  is  a  gain  factor.  As  seen  in  Equa- 
tion (39),  PPF  means  positively  feeding  the  structure  coordi- 
nate back  to  the  compensator  and  positively  feeding  back  the 
compensator  coordinate  times  a  gain  to  the  structure. 

The  system  response  characteristics  of  this  algorithm 
are  shown  in  Figure  9 .  Assuming  that  the  structure  motion  at 
steady  state  is  of  the  form: 


lit)    =  oe 


iut 


(40) 


the   compensator  will    respond  as 


v\(t)    =  pei<wfc-*> 


(41) 


where  the  phase  angle,  <f>,    is  defined  by 


<|>  =  tan 


-i 


2C, 


1  - 


(0 


\2 


(42) 
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Figure  9  Positive  Position  Feedback 


When  the  structure  vibrates  at  a  frequency  much  lower  than  the 
compensator  natural  frequency,  the  phase  angle  approaches 
zero.  Substituting  Equation  (41)  with  0=0  back  into  Equation 
(39a)  gives: 


I   +  2C<o£   +   (co2  -  <?Pg>2)£  =  0 


(43) 


resulting  in  a  decrease  in  the  stiffness  term.  When  the  struc- 
ture and  compensator  have  the  same  natural  frequency,  the 
phase  angle  is  ir/2 .  In  this  case,  the  structural  equation  is 
modified  to: 


I   +  (2C<o  +  <?Pto)  £  +  (o2l   =  0 


(43) 
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showing  an  increase  in  damping.  When  the  structure  is  vibra- 
ting a  frequency  much  greater  than  the  compensator  natural 
frequency,  the  phase  angle  approaches  it.  Making  the  same  sub- 
stitution yields: 

I   +  2C<d£   +    (a)2   +  <?P<o2)£   =  0  (45) 

causing  an  increase  in  the  stiffness  term.  These  equations 
clearly  show  that,  in  order  to  maximize  damping,  coc  must  be  as 
closely  matched  to  u>  as  possible. 

This  analysis  will  go  beyond  the  simple  scalar  case 
and  investigate  a  multi -modal  vibration  suppression  system. 
Specifically,  the  first  three  modes  will  be  controlled  using 
two  different  schemes.  First,  three  collocated  sensor/actuator 
pairs  connected  to  three  compensators  tuned  to  the  first  three 
modal  frequencies  will  be  analyzed.  The  second  system  will 
utilize  only  one  collocated  sensor/actuator  pair  connected  to 
the  same  compensators.  Figures  10  and  11  illustrate  these  two 
configurations  with  simple  block  diagrams. 

In  a  multi -modal  case, . the  system  equations  take  the 
form: 

I   +  Di   +  Kl  =  CTGKx\  (46a) 

ii   +  Z?cf|   +  Kcx\  =  KcCl  (46b) 

For  a  system  with  n  structural  modes  and  m  compensators,  £  is 
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Three  Actuators  and  Three  Compensators  Controlling  Three  Modes 
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Figure   10   Three  Actuator   /   Three   Compensator   Control    Scheme 


Single  Actuator  and  Three  Compensators  Controlling  Three  Modes 
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Figure   11   Single  Actuator/Three   Compensator   Control    Scheme 
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a  [n  x  1]  vector  and  77  is  a  [m  x  1]  vector.  D  and  Dc  are 
[n  x  n]  and  [m  x  m]  diagonal  matrices  containing  the  damping 
terms  2^-to,-  and  2£ciuci.  K  and  Kc  are  [n  x  n]  and  [m  x  m]  diag- 
onal matrices  containing  the  squares  of  the  natural  frequen- 
cies co2-  and  to2u.  G  is  a  [m  x  m]  diagonal  matrix  of  compensator 
gains.  C  is  a  fully  populated  [m  x  n]  participation  matrix 
which  determines  the  influence  of  each  sensor/actuator  pair  on 
each  compensator  and  vice  versa. 

The  formulation  of  these  matrices  is  straightforward 
except  for  the  participation  matrix.  Equating  the  right  hand 
sides  of  Equations  (33)  and  (46a) : 

&TFu   =  CTGKr\  (45) 

and  assuming  the  control  vector  u  is  equal  to  the  vector  GK17, 
C  is  equal  to  (<£TF)T.  u  is  a  [p  x  1]  vector,  where  p  is  the 
number  of  actuators,  which  forces  the  dimensions  of  3>TF  to  be 
[n  x  p] .  A  modification  must  be  made  in  the  case  of  the  single 
actuator  (u  =  1) .  In  order  to  keep  the  G  and  K  matrices  in  the 
same  form  and  ensure  dimensional  compatibility,  a  [1  x  m] 
vector  A  whose  elements  are  unity  must  be  included  such  that: 

<&TFu  =  CTAGKx\  (4  8) 
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Now  the  entire  system  can  be  put  into  the  matrix  form: 


i 

+ 

rD    0 

e 

+ 

A 

0   Dc 

A 

CT(A)GK 


-KC(AT)C 


(49) 


As  was  done  in  Chapter  II.C.l  with  Equation  (38),  a 
state- space  transformation  is  made  to  perform  the  simulations. 
The  same  initial  conditions  are  used  to  excite  the  tri -modal 
vibration  previously  shown  in  Figure  4.  The  compensator  natur- 
al frequencies  were  chosen  to  exactly  match  the  structural 
frequencies  in  order  to  maximize  the  damping  effect.  In  this 
case,  the  com-pensator  damping  ratio,  fc,  should  be  as  small 
as  possible.  As  "shown  in  Figure  12,  this  creates  a  resonance 
peak  near  the  compensator  frequency  allowing  the  largest  gain 
possible  to  act  on  the  affected  mode  in  the  active  damping  re- 
gion. This  will  also  diminish  the  effects  of  active  flexibil- 
ity on  lower  modes.  For  these  simulations,  it  was  assumed  that 
each  compensator  could  be  designed  to  have  a  damping  ratio 
equal  to  . 1.  In  actuality,  it  will  be  very  unlikely  that  the 
compensators  can  be  designed  as  precisely  as  they  have  been 
here  and  a  decrease  in  system  performance  is  expected.  To 
prevent  a  significant  decrease,  the  compensator  damping  ratio 
should  be  increased  so  that  the  active  damping  region  is  wid- 
ened as  much  as  possible. 

Figure  13  shows  the  response  of  the  beam  under  the  in- 
fluence of  three  actuators  at  the  base,  in  the  middle  and  at 
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Figure  12   PPF  Compensator  Transfer  Function  Magnitude 

the  tip.  The  effectiveness  of  this  system  is  seen  in  how 
quickly  all  three  modes  are  suppressed.  A  single  actuator  at 
the  base  also  provides  a  sufficient  degree  of  control  as  seen 
in  Figure  14.  The  effect  of  actuator  placement  is  evident  in 
a  comparison  with  Figures  15  and  16  which  correspond  to  a 
single  actuator  in  the  middle  and  at  the  tip,  respectively, 
feeding  the  same  compensators  with  the  same  gains.  However, 
the  performance  of  these  last  two  configurations  can  be  im- 
proved by  increasing  the  gains.  Due  to  the  fact  that  a  larger 
moment  can  be  placed  at  the  tip  than  at  the  base  without  dri- 
ving the  beam  unstable,  the  allowable  gains  are  two  orders  of 
magnitude  greater  for  Case  3  than  Case  1.  The  individual  modal 
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responses  can  also  be  tailored  by  varying  the  gains  such  that 
the  damping  of  one  mode  is  maximized  in  relation  to  the 
others . 
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Figure  13  PPF  Three  Actuators 

In  the  scalar  case,  system  stability  is  clearly  de- 
fined. A  simple  manipulation  of  the  equations  of  motion  and 
application  of  the  Routh-Hurwitz  criterion  results  in  the 
condition: 


0  <  G  <   1 


(50) 


This  method  is  not  applicable  to  the  much  more  complex  multi- 
modal case.  Fanson  and  Caughey  have  developed  the  following 
stability  condition: 
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Figure  14  PPF  /  Actuator  at  base 
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Figure  15  PPF  /  Actuator  in  middle 
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Figure  16  PPF  /  Actuator  at  tip 


K  -  CTGKC  >   0 


(51) 


for  employing  PPF  in  a  multi -modal  system  [Ref .  8] .  This  con- 
dition is  limited  in  that  this  single  criterion  is  a  function 
of  m  variables  which  are  the  gains  of  the  individual 
compensators.  The  procedure  used  for  these  simulations  was  to 
fix  two  gains  while  varying  the  third  to  determine  the  range 
where  the  eigenvalues  of  the  assembled  matrix  were  all 
positive . 

3.   Strain  Rate  Feedback  (SRF) 

The  final  algorithm  investigated  is  also  defined  in 
modal  form  as: 
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£   +  2Co>£   +  fa)2£   =  -G(x>2i) 


il   +  2Cco)ctl   +  (*2cr\   =  fa)c£ 


(52a) 
(52b) 


where  the  variables  are  the  same  as  those  previously  defined 
for  PPF.  In  this  case,  the  structural  velocity  coordinate  is 
fed  back  to  the  compensator  and  the  compensator  position 
coordinate  is  fed  back  times  a  negative  gain  to  the  structure. 
The  gain  must  be  negative,  as  will  be  shown  in  the  following 
analysis,  so  that  the  system  performance  will  be  optimized 
over  the  most  useful  frequency  range.  Figure  17  graphically 
summarizes  the  performance  characteristics  of  SRF. 

If  we  again  assume  that  the  structural  modal  coordi- 
nate is  of  the  form: 


£(fc)  =  <xeiwt 


(53) 


the  compensator  response  will  be: 


Ti(t)  =  pe^"*-*' 


(54) 


where  <f>   is   now  defined  by: 


4>  =  tan  1 


1  (       > 

0) 

V       c 


-   1 


2C^ 


fa) 


a      J 


(55) 


35 


7t/2-  - 


Til  2 


Active  Postive  Damping 

Active 
Stiffness 

4 ► 

1 

Active  Negative  Damping 

1 

p 

Wr 


PHASE  ANGLE  VS  FREQUENCY 
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Figure  17  Strain  Rate  Feedback 

When  <j>  «  wc/  the  phase  angle,  <f>,  approaches  -it/2.  Substitu- 
ting Equation  (54)  with  <f>  =  -ir/2  back  into  Equation  (52a) 
gives : 


I    +    (2C<o    +  (?Pco)  £    +  G>2£    =  0 


(56) 


resulting  in  an  increase  in  damping.  When  w  =  wc,  <f>   =  0  and  a 
similar  substitution  yields: 


I   +  2C«£  +  (w2  +  <?Po2)£  =  0 


(57) 


increasing  the  stiffness  of  the  structure.  Finally,  when  u>  >> 
wc,  <f>   =  7r/2  and  the  structure  response  is: 
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I   +    (2Cco  -  G$o)£   +  g>2£   =  0 


(57) 


showing  a  decrease  in  damping.  The  obvious  conclusion  is  to 
design  the  compensator  such  that  the  structural  frequencies  of 
the  modes  to  be  damped  fall  well  below  the  compensator  natural 
frequency. 

The  SRF  simulations  were  done  using  the  same  format  as 
the  PPF  simulations.  The  first  three  structural  modes  were 
controlled  using  three  compensators  and  the  same  four  sensor/ 
actuator  configurations.  Additionally,  the  possibility  of  con- 
trolling the  first  three  modes  using  only  one  compensator  and 
one  actuator  is  investigated.  The  participation  matrix  formu- 
lation is  identical  to  the  PPF  case.  The  state-space  transfor- 
mation is  applied  to  the  system  defined  by: 


-a>c(Ar)C  Dc 


K     CT(A)GK 


=   0 


(59) 


Selecting  a  precise  compensator  frequency  is  not  as 
clearly  defined  as  it  is  with  PPF.  Additionally,  the  stability 
condition  is  not  as  clearly  as  defined  as  in  the  PPF  case  due 
to  the  fact  that  the  closed- loop  damping  and  stiffness  matri- 
ces can  not  be  symmetrized.  However,  the  wider  active  damping 
region  gives  the  designer  a  lot  more  flexibility  in  the  cir- 
cuit design.  As  long  as  the  compensator  frequency  is  greater 
than  the  structural  frequency,  a  certain  amount  of  active 
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damping  will  be  provided.  Although  the  stabilizing  feedback 
gain  matrix  is  heavily  dependent  upon  the  structure  and  com- 
pensator frequencies,  some  numerical  analyses  predict  that 
higher  compensator  frequencies  produce  larger  damping  ratios. 
A  drawback  to  frequencies  in  this  range  is  that,  as  shown  in 
Figure  18,  the  magnitude  of  the  compensator  transfer  function 
approaches  zero  limiting  the  amount  of  control.  Further  study 
is  required  on  this  point.  Again,  the  compensators  were  assum- 
ed to  have  damping  ratios  equal  to  .1.  To  highlight  the  ef- 
fects of  varying  compensator  frequencies,  simulations  were  run 
for  the  three  actuator/three  compensator  configuration  with 
two  different  sets  of  wc  at  the  two  extremes.  First,  the  com- 
pensators were  assigned  frequencies  ten  per  cent  higher  than 
each  of  the  structural  modes.  This  response  is  seen  in  Figure 
19.  Although  the  higher  modes  are  controlled,  the  first  mode 
is  not  significantly  damped.  The  system  is  operating  too  far 
in  the  active  stiffness  region.  The  dominance  of  the  active 
stiffness  can  be  seen  in  the  higher  frequency  of  the  first 
mode.  The  second  set  of  compensator  frequencies  were  50  times 
higher  than  each  of  the  structural  modes.  This  has  the  effect 
of  placing  all  three  structural  modes  in  the  active  damping 
region  of  all  three  compensators.  This  response  is  shown  in 
Figure  20.  In  this  case,  the  transfer  function  magnitudes  are 
so  small  that  gains  four  orders  of  magnitude  larger  than  the 
previous  case  are  needed.  With  these  large  gains,  the  higher 
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Figure  18  SRF  Compensator  Transfer  Function  Magnitude 
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Figure  19  SRF  wc  =  1.1  x  u 


39 


TIP  MOTION  -  Tin    SRF  U=3  Gl=1.6  62=G3=1 
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Figure  2  0  SRF  wc-  =  5  0  x  co 

modes  are  almost  immediately  suppressed.  The  first  mode,  hav- 
ing the  lowest  relative  gain,  is  not  controlled  as  efficient- 
ly. A  significant  improvement  in  performance  occurs  as  wc 
moves  farther  away  from  w.  The  same  variation  in  the  system 
response  when  the  single  actuator  moves  from  the  base  to  the 
tip  is  seen  as  in  the  PPF  simulations. 

The  control  system  can  be  simplified  by  only  using  one 
compensator  and  one  actuator  at  the  base  as  shown  in  Figure 
21.  The  motivation  for  this  approach  comes  from  the  results  of 
the  three  compensator  case  where  all  three  structural  frequen- 
cies were  below  the  smallest  compensator  frequency.  Theoreti- 
cally, a  compensator  will  provide  active  damping  to  all  modes 
with  a  lesser  frequency.  For  this  simulation,  the  compensator 
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frequency  was  set  at  ten  times  the  third  modal  frequency.  The 
system  matrices  must  be  modified  because  the  compensator  equa- 
tion is  now  in  a  scalar  form.  As  shown  in  Figure  22,  this 
configuration  effectively  controls  all  three  modes  with  maxi- 
mum effect  on  the  higher  modes.  However,  this  method  will  only 
work  for  a  few  densely  spaced  modes.  If  the  frequency  range  is 
expanded  too  far,  the  lower  modes  will  not  receive  a  strong 
enough  control  input. 

This  approach  was  also  taken  using  an  alternative 
method.  The  sign  of  the  gain  applied  to  the  compensator  ouput 
is  changed  such  that  the  structural  modal  equation  is  now: 

I   +  2C&>£   +  cd2£   =  Gto2T\  (60) 

This  has  the  effect  of  swapping  the  regions  of  active  positive 
damping  and  active  negative  damping  and  making  the  region  cen- 
tered around  wc  active  flexibility.  Now,  the  objective  is  to 
set  the  compensator  frequency  below  the  first  mode  placing  all 
the  modes  in  the  active  positive  damping  region.  In  this  simu- 
lation, the  compensator  frequency  was  set  at  80%  of  the  first 
modal  frequency.  A  first  attempt  was  made  using  the  same  com- 
pensator damping  ratio  (fc  =  .1)  .  This  was  an  unstable  system 
due  to  the  resonance  peak  of  the  compensator  transfer  function 
occuring  in  the  active  flexibility  region.  To  improve  the  re- 
sponse, the  damping  ratio  of  the  compensator  had  to  be  in- 
creased such  that  there  was  no  resonance  and  the  flexibility 
effect  was  diminished  with  respect  to  the  damping  effect.  A 
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SINGLE  ACTUATOR  AND  SINGLE  COMPENSATOR  CONTROLLING  THREE  MODES 

SENSOR 

INPUT          OUTPUT 

1 

Compensator  — 

SRF 

ACTUATOR 

Figure  21  Single  Actuator/Single  Compensator  Control  Scheme 


value  for  fc  of  .95  produced  the  response  in  Figure  23.  As  ex- 
pected, the  higher  modes  are  not  damped  out  as  quickly  as  in 
the  previous  case.  To  prevent  the  first  mode  from  being  driven 
unstable,  the  gain  must  be  held  relatively  low  which  reduces 
the  impact  at  higher  frequencies.  In  practice,  the  effective 
bandwidth  of  the  actual  compensator  must  be  considered.  The 
compensator  will  damp  only  those  modes  densely  packed  within 
the  range  of  control. 

Deriving  a  stability  condition  for  SRF  similar  to 
Equation  (51)  for  PPF  proved  difficult.  A  transformation  to 
symmetrize  the  modal  equations  could  not  be  found.  In  their 
unsymmetric  form,  the  equations  could  not  be  manipulated  to 
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Figure  22  SRF  Single  Actuator  /  Negative  Gain 
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Figure  23  SRF  Single  Compensator  /  Positive  Gain 
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form  a  simple  stability  condition.  The  alternate  method  used 
was  to  ensure  all  the  eigenvalues  of  the  system  matrix  of  the 
state- space  equation  were  negative.  Again,  two  gains  were  held 
fixed  while  varying  the  third  to  determine  the  stable  region. 
4.   Comparison  of  PPF  and  SRF 

A  strict  side-by- side  comparison  of  these  two  methods 
is  not  feasible  due  to  the  fact  that  the  compensators  are  not 
performing  the  same  functions  in  both  methods.  Also,  since  the 
stable  gain  regions  were  so  varied,  a  uniform  set  of  operating 
parameters  that  provided  sufficient  control  in  both  systems 
could  not  be  defined. 

Positive  Position  Feedback  is  preferable  for  applica- 
tions where  a  single  mode  needs  to  be  quickly  damped.  Control 
can  be  maximized  by  taking  advantage  of  the  compensator  reso- 
nance occurring  in  the  active  damping  region.  PPF  also  has  the 
advantage  of  using  the  sensors  to  measure  the  structural  dis- 
placement rather  than  the  more  difficult  measurement  of  struc- 
tural velocity  as  is  needed  with  SRF.  For  multi-modal  control, 
the  frequency  requirements  are  too  restrictive.  One  compensa- 
tor per  mode  is  required  in  this  case.  Unless  the  structural 
characteristics  are  uniform  and  precisely  known  and  unless  the 
compensator  is  closely  matched  with  the  modal  frequency,  PPF 
will  be  operating  at  less  than  optimum. 

Strain  Rate  Feedback  has  the  advantage  of  controlling 
modes  over  a  wider  range  of  frequencies.  As  was  demonstrated, 
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a  single  compensator  feeding  a  single  actuator  can  control 
three  modes.  This  is  important  in  terms  of  reducing  circuit 
complexity  and  manufacture  and  implementation  costs.  However, 
the  active  damping  occurs  in  a  region  where  the  compensator 
gain  is  not  at  a  maximum.  This  has  the  effect  of  limiting  the 
amount  of  control  on  modes  far  from  the  compensator  frequency. 

The  robustness  of  these  two  control  schemes  is  demon- 
strated in  Figures  24  and  25.  The  scenario  simulated  here  is 
that  the  actual  structural  frequencies  experienced  on-orbit 
are  10%  less  than  the  predicted  values  around  which  the  com- 
pensators were  designed.  All  other  parameters  are  kept  con- 
stant. This  new  response  of  the  "modified  structure"  is  shown 
with  the  dashed  line  and  is  compared  to  the  response  of  the 
"original  structure"  shown  with  the  solid  line. 

PPF  produces  almost  an  identical  response  for  both 
systems  with  only  miniscule  differences.  Initial  expectations 
were  that  this  frequency  variation  would  produce  a  much  more 
significant  difference.  Further  study  revealed  that  the  ro- 
bustness of  the  system  is  very  dependent  on  the  compensator 
damping  ratio.  Small  values  produce  a  phase  angle  curve  that 
approaces  a  square  wave  which  were  seen  to  emphasize  the  ad- 
verse effects  of  a  frequency  change.  Likewise,  large  values 
stretch  out  the  phase  angle  curve  allowing  a  less  significant 
change  in  the  system  response.  For  this  particular  model,  the 
chosen  value  is  large  enough  to  prevent  PPF's  frequency  depen- 
dency from  degrading  the  system  response.  The  key  to  providing 
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Figure  24  PPF  Structural  Frequencies  Varied  101 
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Figure  25  SRF  Structural  Frequencies  Varied  10% 
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sufficient  robustness  is  to  have  a  moderate  compensator  dam- 
ping ratio.  Small  values  will  limit  the  bandwidth  of  the  com- 
pensator and  large  values,  while  still  being  robust,  will 
decrease  the  magnitude  of  the  damping  provided. 

A  variation  of  only  the  first  mode  is  seen  when  using 
SRF.  The  "modified"  first  modal  frequency  has  moved  closer  to 
the  edge  of  the  bandwidth  of  the  compensators  such  that  the 
amount  of  damping  able  to  be  provided  has  decreased  slightly. 
The  structure  is  also  more  flexible  as  seen  in  the  lower  fre- 
quency. SRF  is  much  less  dependent  on  the  compensator  damping 
ratio  in  terms  of  robustness.  The  system  is  operating  in  a 
region  where  the  slope  of  the  phase  angle  curve  is  not  signif- 
icantly altered  by  variations  of  the  damping  ratio.  Simula- 
tions showed  that  the  system  was  equally  robust  with  damping 
ratios  different  by  two  orders  of  magnitude. 
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III.   EXPERIMENTAL  ANALYSIS 

The  characteristics  of  Positive  Position  Feedback  will  be 
studied  using  the  Naval  Postgraduate  School's  Flexible  Space- 
craft Simulator  (FSS) . 

A.   PHYSICAL  SET-UP 

The  mechanical  lay-out  of  the  experiment  is  shown  in  Fig- 
ure 26.  Reference  4  details  the  circuit  design  and  the  choice 
and  placement  of  the  piezoceramic  sensors  and  actuators  which 
comprise  the  active  damping  control  system.  A  detailed  view  of 
the  sensor  and  actuator  locations  is  shown  in  Figure  27.  The 
ceramics  are  mounted  on  both  sides  of  the  flexible  arm  with 
opposite  polarities  such  that  the  sensor  voltages  and  actuator 
strains  are  consistent.  The  damping  circuit  consists  of  three 
sections:  a  charge  amplifier,  a  low-pass  filter  and  a  power 
amplifier.  The  components  of  the  filter  were  chosen  such  that 
the  compensator  values  are  wc  =  .15  Hz  and  £"c  =  .1. 

For  this  analysis,  the  main  body  of  the  simulator  is  held 
fixed  and  the  flexible  arm  is  the  only  component  affected  by 
external  forces.  The  flexible  arm  is  supported  by  air  pads 
which  allow  two-dimensional  motion  in  an  approximately  fric- 
tionless  environment. 
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Figure  26  Flexible  Spacecraft  Simulator  (FSS 
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Figure  27  Piezoceramic  Locations 
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B .   PROCEDURE 

The  objective  of  the  experiment  was  to  investigate  the 
effectiveness  of  the  PPF  algorithm  over  a  range  of  structural 
frequencies  and  compare  the  results  with  theory. 

To  achieve  a  variety  of  frequencies,  the  configuration  of 
the  flexible  arm  was  changed  by  attaching  and  detaching  the 
concentrated  masses.  The  natural  frequency,  co,  for  a  simple 
spring-mass  system  is  defined  by  the  relationship: 


0)  = 


\ 


k  (61) 

m 


where  m  is  the  mass  and  k  is  the  stiffness.  The  masses  are 
attached  to  the  arm  such  that  the  stiffness  remains  constant. 
Therefore,  adding  mass  will  lower  the  frequency  from  the  ini- 
tial value  and  vice  versa.  The  selection  of  where  to  make  mass 
adjustments  was  a  trial  and  error  procedure  until  frequencies 
of  sufficient  range  and  separation  were  found. 

Data  was  obtained  in  the  following  manner.  The  signal  from 
the  sensors  was  sent  via  a  Gould  5900  Signal  Conditioner  to  a 
Gould  RS3800  Stripchart  Recorder.  This  allowed  obtaining  a 
time  history  of  sufficient  length  needed  to  calculate  damping 
ratios.  Phase  angle  measurements  were  made  using  a  Hewlett- 
Packard  54601A  Oscilloscope  which  displayed  both  the  signal 
from  the  sensors  and  the  signal  to  the  actuators.  Frequency 
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measurements  were  made  using  both  sets  of  equipment.  The  data 
was  averaged  over  several  run  for  each  configuration. 

With  only  a  single  compensator  in  the  feedback  circuit, 
the  system  will  only  be  effective  on  the  mode  for  which  it  is 
designed.  Therefore,  in  this  experiment,  only  the  first  mode 
of  the  flexible  arm  was  excited  so  that  the  cleanest  and 
smoothest  signal  could  be  recorded. 

The  damping  ratio,  £ ,  was  calculated  using  the  log  decre- 
ment method: 


2nn         \Afj 


(62) 


where  Aj  is  the  initial  amplitude,  Af  is  the  final  amplitude 
and  n  is  the  number  of  cycles  between  those  two  amplitudes. 
The  natural  frequency  is  computed  using  the  damping  ratio  and 
the  damped  period,  rd: 


_2n 

'n 


<■>-  =  " (63) 


x. 


\/T^T 


However,  the  damping  ratios  are  small  enough  with  respect  to 
the  frequency  measurement  accuracy  that  we  can  assume  the  nat- 
ural frequency  to  be  the  same  as  the  damped  frequency.  The 
phase  angle  was  calculated  from  the  relationship: 

<t>  =  g>  At  (360°/cycJe)  (64) 
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where  w  is  the  structural  frequency  in  Hz  and  At  is  the  time 
difference  between  the  sensor  and  actuator  signals  as  seen  on 
the  oscilloscope.  The  theoretical  value  for  the  phase  angle  at 
each  structural  frequency  is  calculated  using  Equation  42. 

C .   RESULTS 

As  shown  in  the  plot  in  Figure  28,  data  was  obtained  in 
all  three  regions  of  active  flexibility,  damping  and  stiff- 
ness. (Refer  to  Figure  9).  Although  the  experimental  values 
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Figure  2  8  FSS  Flexible  Arm  Phase  Angles 

for  phase  angle  vary  from  theory  (see  Table  I)  as  much  as 
seven  degrees  in  the  worst  case  and  considering  the  +/-  2 
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degree  measurement  accuracy,  the  data  clearly  follows  the 
trend  we  expect  and  we  can  assume  the  circuit  is  functioning 
properly. 
Table  I 


FLEXIBLE  SPACECRAFT  SIMULATOR  -  PPF 

CASE 

PHASE  ANGLE  (Degrees) 

Theoretical 

Experimental 

1 

19 

15 

2 

29 

29 

3 

32 

36 

4 

55 

59 

5 

90 

84 

6 

111 

107 

7 

139 

146 

8 

150 

156 

The  parameters  describing  the  dynamic  response  of  the 
flexible  arm  are  listed  in  Table  II.  Examining  the  change  in 
frequency  from  the  freely  vibrating  (open- loop)  to  the  active- 
ly damped  (closed-loop)  cases,  the  expected  effect  on  the 
stiffness  term  is  evident.  At  the  lower  frequencies,  a  de- 
crease is  observed  indicating  the  arm  has  become  more  flexi- 
ble. At  and  near  the  compensator  frequency,  there  is  little  or 
no  change.  The  higher  frequencies  exhibit  an  increase  meaning 
the  arm  has  become  stiffer. 
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Table  II 


FLEXIBLE  SPACECRAFT  SIMULATOR  -  PPF 

CASE 

FREQUENCY  (Hz) 

DAMPING  RATIO 

Free 

Damped 

Free 

Damped 

%  A 

1 

.116 

.113 

.0095 

.0070 

-27 

2 

.128 

.125 

.0104 

.0114 

10 

3 

.133 

.128 

.0077 

.0243 

215 

4 

.141 

.140 

.0059 

.0325 

450 

5 

.150 

.150 

.0050 

.0221 

345 

6 

.154 

.156 

.0040 

.0158 

299 

7 

.165 

.168 

.0056 

.0145 

161 

8 

.177 

.180 

.0035 

.0093 

165 

A  very  pronounced  variation  is  seen  in  the  change  in  the 
damping  ratio  when  the  feedback  system  is  activated.  At  the 
lowest  observed  frequency,  the  arm  has  become  so  flexible  that 
the  damping  worsens  by  27  per  cent.  There  is  a  steady  improve- 
ment as  the  structure  moves  into  the  active  damping  region. 
The  maximum  effect  is  observed  with  Case  4  showing  a  45  0  per 
cent  improvement .  The  magnitude  of  the  control  provided  is 
best  appreciated  visually  as  seen  in  Figure  29.  Even  at  the 
higher  frequencies,  while  not  in  the  active  damping  region,  an 
improvement  in  the  damping  ratio  is  produced  by  the  increased 
stiffness  of  the  arm. 
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One  discrepancy  evident  in  this  data  is  that  the  greatest 
improvement  in  the  damping  provided  by  the  compensator  does 
not  occur  at  a  structural  frequency  equal  to  the  compensator 
frequency.  The  cause  of  this  variation  has  been  difficult  to 
pin-point.  In  order  for  the  transfer  function  resonance  peak 
to  occur  at  .141  Hz,  where  the  maximum  increase  was  observed, 
the  compensator  damping  ratio  would  have  to  increase  to  .34  or 
the  compensator  frequency  would  have  to  decrease  to  .142  Hz. 
However,  this  is  not  consistent  with  the  phase  angle  data  pre- 
viously mentioned.  This  rules  out  any  small  variations  in  the 
feedback  circuit  as  possible  sources  of  error.  There  are  some 
physical  factors  which  may  have  caused  inconsistencies  between 
the  eight  configurations.  Other  experiments  were  being  conduc- 
ted while  this  one  was  in  progress  which  did  not  allow  a 
uniform  orientation  of  the  flexible  arm  for  all  of  the  config- 
urations. Varying  frictional  or  gravitational  effects  may  have 
adversely  affected  the  consistency  of  the  results.  Further  de- 
tailed experimental  work  is  needed  to  precisely  determine 
where  the  maximum  effect  actually  occurs. 

This  data  highlights  the  advantages  and  disadvantages  of 
the  Positive  Position  Feedback  algorithm.  When  the  compensator 
is  precisely  tuned  to  the  structural  frequency,  the  system  is 
operating  near  the  maximum  gain  of  the  compensator  transfer 
function.  This  enables  the  active  damping  control  to  have  its 
greatest  effect.  However,  any  deviation  from  this  optimum  re- 
gion reduces  the  impact  of  the  control  system.  As  seen  in  this 
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Figure  29  Flexible  Arm  First  Mode  Response 
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experiment,  only  a  25  per  cent  reduction  in  the  structural 
frequency  results  in  the  compensator  degrading  the  dynamic  re- 
sponse of  the  system.  If  the  frequencies  can  not  be  matched 
precisely,  it  is  better  to  err  into  the  active  stiffness  re- 
gion where  the  structural  characteristics  are  not  adversely 
affected.  The  efficiency  of  the  control  law  with  a  bounded 
gain  space,  0<G<1,  is  heavily  dependent  upon  the  compensator 
parameters  coc  and  fc.  It  is  recommended  that  moderately  higher 
compensator  damping  ratios  be  used  so  that  a  relatively  large 
amount  of  damping  is  provided  while  still  maintaining  suffi- 
cient robustness. 
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IV.   CONCLUSIONS 

Piezoceramic  sensors  and  actuators  have  been  shown  to  pro- 
vide an  effective  means  of  controlling  structural  vibrations. 
An  increase  of  450  per  cent  in  the  damping  ratio  of  the  flexi- 
ble arm  of  the  Flexible  Spacecraft  Simulator  was  demonstrated 
using  the  Positive  Position  Feedback  algorithm.  This  system 
provides  the  greatest  degree  of  control  when  the  compensator 
is  precisely  tuned  to  the  modal  frequency.  Its  main  limitation 
is  that  only  a  single  mode  can  be  controlled  per  compensator. 

Positive  Position  Feedback  and  Strain  Rate  Feedback  have 
significantly  different  control  characteristics  which  were 
analytically  demonstrated  for  a  cantilevered  beam.  SRF  is  more 
advantageous  than  PPF  in  the  multi-modal  environment.  With  a 
wider  active  damping  region,  it  has  the  ability  to  stabilize 
more  than  one  mode  given  a  sufficient  bandwidth. 

Recommendations  for  further  study  include  a  number  of 
topics.  First,  the  dynamic  response  of  the  FSS  flexible  arm 
could  be  improved  by  incorporating  additional  compensators 
into  the  current  configuration.  The  first  mode  is  very  rarely 
the  only  mode  excited  during  slew  maneuvers  and  the  higher 
modes  should  be  considered.  This  problem  may  be  more  easily 
solved  by  building  a  Strain  Rate  Feedback  circuit  and  integra- 
ting it  with  the  sensors  and  actuators  currently  mounted  on 
the  arm.  Also,  the  effect  of  multiple  sensors  and  actuators  on 
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the  arm  can  be  investigated  with  an  optimization  study  of  sen- 
sor/actuator placement.  There  are  a  variety  of  control  laws 
which  can  be  used  in  this  application.  As  an  example,  a  combi- 
nation of  PPF  and  SRF  could  be  developed  which  may  incorporate 
the  strengths  of  each  and  cancel  out  the  weaknesses. 

The  Flexible  Spacecraft  Simulator  is  an  excellent  tool  for 
studying,  not  only  active  damping  control  techniques,  but 
other  spacecraft  attitude  control  and  dynamics  problems.  Cur- 
rent research  endeavors  include  designing  an  improved  thruster 
for  attitude  control  and  incorporating  a  robotic  arm  onto  the 
main  body. 
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APPENDIX 

%  Simulation  of  Beam  Vibration  using  Finite  Element  Method 

%  Transformation  into  modal  form  using  PPF  with  three  control 

laws 

%  utilizing  three  actuators  on  elements  1,5,8. 

% 

%  Define  physical  constants: 


n=8;  %  number  of  elements 

h=.l;  %  element  length 

d= . 5 ;  %  element  density 

EI=.9;  %  modulus  of  elasticity  times  moment  of 

inertia 

% 

%  Element  mass  matrix: 

% 

Mul=(d*h/420) *  [156    22*h;    22*h   4*h^2]  ; 

Mur= (d*h/420) *  [54    -13*h;    13*h    -3*h*2] ; 

Mll=Mur' ; 

Mlr=(d*h/420) *  [156    -22*h;     -22*h   4*hA2]  ; 

o, 

"S 

%  Element  stiffness  matrix: 

% 

Kul=(EI/hA3) * [12  6*h;  6*h  4*h^2] ; 

Kur=(EI/hx3) * [-12  6*h;  -G*h  2*hA2] ; 

Kll=Kur' ; 

Klr= (EI/h^3) * [12  -6*h;  -6*h4*h^2]; 

% 

%  Combine  matrices  for  n  element  system.  Disregard  first  two 

rows 

%  and  columns  for  clamped  end. 

% 

for  i=l:n-l, 

M(2*i-l:2*i,2*i-l:2*i)=Mul+Mlr; 

K(2*i-l:2*i,2*i-l:2*i)=Kul+Klr; 
% 

M(2*i-l:2*i,2*i+l:2*i+2)=Mur; 

K(2*i-l:2*i,2*i+l:2*i+2)=Kur; 


M(2*i+l:2*i+2,2*i-l:2*i)=Mll; 
K(2*i+l:2*i+2,2*i-l:2*i)=Kll; 


end 


M(2*n-l:2*n,2*n-l:2*n)=Mlr; 

K(2*n-l:2*n,2*n-l:2*n) =Klr; 

% 

Mi=inv(M);  %  Inverse  of  mass  matrix. 
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%  Obtain  eigenvalues  and  eigenvectors: 

o 

[phi, lambda] =eig (K,M) ; 
% 

%  Extract  natural  freqencies: 
% 

w=diag (lambda) ; 
% 

%  Sort  from  low  to  high  frequency: 
% 

[wl , num] =sort (w) ; 
w=w  (num)  ; 
for  j=l:2*n 

lambda  (j  ,  j)  =w(j  )  ; 
end 
% 

%  Sort  eigenvectors: 
% 

phi=phi ( : , num) ; 
% 

%  Calculate  generalized  masses: 
% 

Mg=phi ' *M*phi ; 
% 

%  Calculate  weighted  modal  matrix: 
% 
for  i=l:2*n 

phiw ( : , i ) =phi ( : , i ) /sqrt (Mg ( i , i )  )  ; 
end 

Mw=phiw' *M*phiw; 
Kw=phiw' *K*phiw; 
% 

%  Actuator  location  matrix  (actuators  on  elements  1,5,8): 
% 

F=zeros (16, 3) ; 

F (2,1) =1;F (8,2) =-l;F (10,2) =1;F( 14,3) =-l;F( 16,3) =1; 
Ft=phiw' *F; 
% 

%  Assign  structural  modal  damping  values: 
% 

z (1)  =  . 0001; z  (2) =.0004; z (3) =.0009; z (4)  =  .0015; z (5)  =  .0037; z (6) - 
.0089; 

z(7)=.01;z (8) =.025;z(9) =.042;z  (10) = . 067; z (11) = . 098  ;  z (12)  =  .  11; 
z(13)=.35;z(14)=.86;z(15)=1.21;z(16)=1.77; 
% 

%  Assemble  structural  damping  matrix: 
% 

D=zeros (2*n) ; 
for  i=l:2*n 

D(i,i)=2*z(i)*sqrt (Kw(i,i) ) ; 
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end 
% 

%  Truncate  to  first  three  modes: 
% 

Kwl=Kw(l:3, 1:3) ; 
Ftl=Ft (1:3,1:3) ; 
D1=D(1:3,1:3) ; 
% 

%  Assign  compensator  values: 
% 

wcl=sqrt (Kwl (1,1) ) ;wc2=sqrt (Kwl (2,2) ) ;wc3=sqrt (Kwl (3,3) ) ; 
zcl=. 1; zc2= . 1; zc3= . 1; 
G1=.2  03  8;G2=.01;G3=.01; 
Kc= [wclA2  0  0;  0  wc2^2  0;  0  0  wc3^2]  ; 
Dc=  [2*zcl*wcl  0  0;  0  2*zc2*wc2  0;  0  0  2*zc3*wc3] ; 
G=[G1  0  0;  0  G2  0;  0  0  G3]  ; 
% 

%  Combine  structure  and  compensator  equations  into  matrix 
form: 
% 

Ds= [Dl  zeros (3);  zeros (3)  Dc]  ; 
Ks=[Kwl  -Ftl*G*Kwl;  -Kc*Ftl'  Kc] ; 
% 

%  Put  closed-loop  system  in  state  space  form: 
% 

A=  [zeros (6)  eye (6);  -Ks  -Ds] ; 
B=zeros (12, 1) ; 
C=eye  (12)  ; 
D=zeros (12 , 1) ; 
T-(0: .01:9.99) ; 
U=zeros (T) ; 
X0=zeros (2*n, 1) ; 
% 

%  Initial  displacements  to  excite  3  modes 
% 

X0 ( 1 , 1 )  =  . 00  8 ; X0 (2 , 1 ) = . 05 ; X0 ( 3 , 1 ) = . 01 ; X0 (4 , 1 ) =0 ; X0 ( 5 , 1 ) = . 004 ; 
X0(6,l)=- .086;X0 (7,1)=- . 007 ;X0 (8 , 1) =- . 086 ;X0 (9 , 1) =- . 016 ;X0 (1 
0,1)=0; 

X0(11,1)=- .015;X0(12,1)=.01;X0 (13 , 1) =- . 001 ;X0 (14, 1)=. 173; 
X0(15,1)=.02;X0 (16 , 1) = . 173 ; 

Z0=inv(phiw) *X0;  %  Initial  displacement  in  modal 

coordinates 
Z01  =  Z0  (1:3)  ; 

M0=  [Z01'     00000000    0]'; 
[Y,X] =lsim(A,B,C,D,U,T,M0) ; 

o, 
o 

%  Put  open- loop  system  in  state- space  form: 

% 

Aol= [zeros (3)  eye (3);  -Kwl  -Dl] ; 

Bol=zeros (6,1) ; 

Col=eye (6) ; 
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Dol=zeros (6,1) ; 
MOol  = [Z01'  OOO]; 

[Yol,Xol] =lsim(Aol, Bol , Col ,Dol , U, T,MOol) ; 
% 
% 

%  Convert  back  to  tip  physical  coordinate: 
% 

phiq=phiw (1:16,1:3) ; 
for  i=l:1000 

tm=[Y(i,l)  Y(i,2)  Y(i,3)] ' ; 

q=phiq*tm; 

tip(i)=q(15) ; 

tmol=[Yol (i,l)  Yol(i,2)  Yol(i,3)]'; 

qol=phiq*tmol ; 

tipol (i) =qol (15) ; 
end 
% 
%  Plot  data 

o, 
"6 

plot  (T,tip(:)) 

xlabel  ('TIME,  sec') 

ylabel  ('  TIP  DISPLACEMENT'  ) 

title  ('TIP  MOTION  -  FEM  PPF  U=3  Gl=.2038  G2=G3=.01' 

pause 

plot  (T,tipol  (:)") 

xlabel  ('TIME,  sec') 

ylabel  ( 'TIP  DISPLACEMENT' ) 

title  ('TIP  MOTION  -  FEM  OPEN  LOOP') 

pause 
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